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Abstract 

We review some contributions on fusion rules that were inspired by the work of Sharp, 
in particular, the generating-function method for tensor-product coefficients that he de- 
veloped with Patera. We also review the Kac- Walton formula, the concepts of threshold 
level, fusion elementary couplings, fusion generating functions and fusion bases. We try to 
keep the presentation elementary and exemplify each concept with the simple sti(2)fc case. 
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1. Introduction 



The Patera-Sharp generating function for su(2) tensor products reads [1]: 

^\L,M,N) = {l _ LKm _\ N){l _ MN) (1.1) 

where the multiplicity of the representation (n) in the tensor product (£) <8> (m) is the 
coefficient of L^M m iV n in the series expansion of (1.1). The derivation of this expression 
is based on manipulations of the character generating functions [2] . 

Key concepts on fusion rules have been obtained by looking for the affine-fusion ex- 
tension of this simple-looking expression and its simplest higher-rank relatives. These are: 
the threshold level, fusion elementary couplings and fusion bases. Before reviewing these 
results, we briefly discuss the Kac- Walton formula. This last result can also be linked, 
albeit loosely, to Bob Sharp. Indeed, it is an affine extension of the Racah-Speiser al- 
gorithm for computing tensor-product coefficients, one of Sharp's favorite techniques. He 
presented it in his course on group theory, where two of the authors (PM and MW) learned 
the fundamentals of this subject. 



2. Fusion rules: the set up 

Fusion rules give the number of independent couplings between three given primary 
fields in conformal field theories. We are interested in those conformal field theories having 
a Lie algebra symmetry. These are the Wess-Zumino-Witten models [3] , whose spectrum- 
generating algebra is an affine Lie algebra at integer level. Their primary fields are in 1-1 
correspondence with the integrable representations of the appropriate affine Lie algebra at 
level k. Denote this set by P.! and a primary field by the corresponding affine weight A. 
Fusion coefficients v are defined by the fusion product 

To simplify the presentation, we consider only the algebra su(N). 
An affine weight may be written as 

N-l 

A = ^ XiOJi = [A , Ai, Ajv-i] (2.2) 

i=0 



where Qi denote the fundamental weights of su (AT). If the Dynkin labels \{ are nonnegative 
integers, the weight A is the highest weight of an integrable representation of su(N) at level 
/c, with k defined by k = X^o 1 To tne a ffi ne weight A, we associate a finite weight A 
of the finite algebra su(N): 

N-l 

A = XjiOj = (Ai, Ajy-i) (2.3) 
i=i 

where are the fundamental weights of su(N). Thus A is uniquely fixed by A and k. 



3. Fusion rules and tensor products: the Kac- Walton formula 

The fusion coefficient AT™ v is fixed to a large extent by the tensor-product coefficient 
pertaining to the product of the corresponding finite representations. We denote by 
the multiplicity of the representation v in the tensor product A ® p: 

By abuse of notation, we use the same symbol for the highest weight and the highest- 
weight representation. P+ represents the set of integrable finite weights. The precise 
relation between tensor-product and fusion-rule coefficients is given by the Kac- Walton 
formula [4,5]: 

= £ ^e(w) (3.2) 

■wew , w-i=oeP^ 

w is an element of the affine Weyl group W, of sign e(w), and the dot indicates the shifted 
action: w ■ A = w(X + p) — p, where p stands for the affine Weyl vector p = Xlilo 1 

The Kac- Walton formula can be transcribed into a simple algorithm: one first calcu- 
lates the tensor product of the corresponding finite weights and then extends every weight 
to its affine version at the appropriate level k. Weights with negative zeroth Dynkin label 
are then shift-reflected to the integrable affine sector. Weights that cannot be shift-reflected 
to the integrable sector are ignored (this is the case, for example, for those with zeroth 
Dynkin label equal to —1). 



Here is a simple example: consider the su(2) tensor-product 

(2) ® (4) = (2) © (4) © (6) (3.3) 

and its affine extension at level 4: 

[2, 2] x [0,4] = [2, 2] + [0,4] + [-2, 6] (3.4) 

The last weight must be reflected since it is not integrable: the shifted action of sq, the 
reflection with respect to the zeroth affine root, is 

s ■ [-2, 6] = s ([-2, 6] + [1, 1]) - [1, 1] = [0, 4] (3.5) 

and this contributes with a minus sign (e(so) = — l)j cancelling then the other [0, 4] repre- 
sentation; we thus find: 

[2,2] x [0,4] = [2,2] (3.6) 

The relation between tensor products and fusion was further explored in [6], published 
in a special volume of the Canadian Journal of Physics dedicated to Prof. R.T. Sharp. 

4. The idea of threshold level 

The result of the above computation is manifestly level dependent. Let us reconsider 
the same product, but at level 5. The affine extension of the product becomes 

[3, 2] x [1,4] = [3, 2] + [1,4] + [-1,6] (4.1) 

The last weight is thus ignored and the final result is [3, 2] x [1, 4] = [3, 2] + [1, 4]. For k > 5 
it is clear that there are no truncations, hence no difference between the fusion coefficients 
and the tensor products. Moreover, we see that the representation (4) occurs at level 5 
and higher. We then say that its threshold level, denoted by k , is 5. The threshold level 
is thus the smallest value of k such that the fusion coefficient v is non-zero, when 

Apt 

G {0, 1}, for all levels fcJil If we indicate the threshold level by a subscript, we can 

write 

(2)©(4) = (2) 4 ©(4) 5 ©(6) 6 (4.2) 

1 More generally, we say there are Mx^ couplings, each having its own threshold level ko- 
For fixed {A, fi; is} then, one gets a multi-set of threshold levels. A simple example: su(3) with 
A = fx = v labelling the adjoint representation has threshold levels {2, 3}. 



To read off a fusion at fixed level k, we only keep terms with index not greater than k. 
This implies directly the inequality 



which in turn yields 



AA. (fc) < A/1 fc+1) (4.3) 

Xp, Xp, 



lim N^ P =Nx^ (4.4) 



The concept of threshold level was first introduced in [7] . Its origin is directly rooted 
in the generating-function method applied to fusion rules. This is reviewed in the next 
section, focusing again on the simple su(2) case. 



5. Fusion generating functions 

The result (1.1) on the su(2) tensor-product generating function can be understood 
as follows: all couplings can be described by appropriate products of three tensor-product 
elementary couplings: 

E x =LM E 2 = LN E 3 = MN (5.1) 
qsu{2) can ftms be written compactly as 

G su{2 \L,M,N) = f[—l— (5.2) 

i=l ^ % > 

How can we construct the afhne extension of this generating function? One certainly 
needs to introduce a further dummy variable, say d, in order to keep track of the extra 
variable k. Then one could try to introduce factors of d appropriately. A natural guess is 

(1 — — dLM)(l — dLN)(l — dMN) 1 ' } 

That turns out to be the right answer: this reproduces the su(2)k fusion rules. The 
prefactor is justified as follows: the fusion of the 'vacuum' with itself, [k, 0] x [k, 0] = [k, 0], 
exists at every level and this is precisely taken into account by the factor 1/(1 — d). 
With the threshold level insight, we can also naturally justify the factors of d multiplying 
the three elementary couplings: their power yields their threshold level. This expression 



was first proved in [7]. From this example and that of the su(3) generating function, 
we conjectured that any fusion coupling is characterized by a threshold level. This was 
subsequently checked with the sb(5) case [8]. This is now understood to be a consequence 
of the Gepner-Witten depth rule [3] , as shown in [9] . 

From the above expression, we also infer the existence of fusion elementary couplings. 
For su(2), there are four of them 

E :d: (0) <g> (0) D (0)i E t : dLM : (1) <g> (1) D (0)i, 

(5-4) 

E 2 : dLN : (1) <g> (0) D (l)i, E 3 : dMN : (0) <g> (1) D (l)i. 

As explained above, subscripts indicate the threshold level. 

A re-derivation of (5.3) was presented in [10]. The method used there was amenable to 
generalization, unlike the original proof in [7]. Consequently, [10] displays further examples 
of fusion generating functions. 



6. Tensor products, linear inequalities and elementary couplings 

There are simple combinatorial methods that can be used for calculating su(N) tensor 
products, for instance, the Littlewood-Richardson (LR) rule. It is thus natural to ask 
whether we can read off the threshold level of a coupling from its LR tableau. 

Integrable weights in su(N) can be represented by tableaux: the weight 
(Ai, A2, • • • , Aat_i) is associated to a left justified tableau of N — 1 rows with Ai + A2 + . . . + 
Xn-i boxes in the first row, A2 + . . . + Aat_i boxes in the second row, etc. Equivalently, 
the tableau has Ai columns of 1 box, A2 columns of 2 boxes, etc. The scalar representation 
has no boxes, or equivalently, any number of columns of iV boxes. For instance: 



su(3):(l,l) 



su(4) : (2,3,0) 



(6.1) 



The Littlewood-Richardson rule is a simple combinatorial algorithm that calculates 
the decomposition of the tensor product of two su(N) representations A <8> \i. The second 
tableau (//) is filled with numbers as follows: the first row with l's, the second row with 
2's, etc. All the boxes with a 1 are then added to the first tableau according to following 
restrictions: (1) the resulting tableau must be regular: the number of boxes in a given row 



must be smaller or equal to the number of boxes in the row just above; (2) the resulting 
tableau must not contain two boxes marked by 1 in the same column. All the boxes marked 
by a 2 are then added to the resulting tableaux according to the above two rules (with 
1 replaced by 2) and the further restriction: (3) in counting from right to left and top 
to bottom, the number of l's must always be greater or equal to the number of 2's. The 
process is repeated with the boxes marked by a 3, 4, ... , N — l, with the additional rule that 
the number of i's must always be greater or equal to the number of i + l's when counted 
from right to left and top to bottom. The resulting Littlewood-Richardson (LR) tableaux 
are the Young tableaux of the irreducible representations occurring in the decomposition. 

Here is a simple su(3) example: (1, 1) ® (1, 1) D 2(1, 1) since we can draw two LR 
tableaux with shape (1,1) and an extra column of three boxes (the total number of boxes 
being preserved, the resulting LR tableau must have 6 boxes): 







1 








1 


2 




2 


1 


(6.2) 



These rules can be rephrased in an algebraic way as follows [11]. Define n^- to be the 
number of boxes i that appear in the LR tableau in the row j. The LR conditions read: 

fc-l k 

Aj-i + n i3-i ~ Yl n v ^ l<k<j<N j^l (6.3) 

i=i i=i 

and 

k k 

- J^ny > 2 < i < k < N and i < N — 1. (6.4) 

The weight \i of the second tableau and the weight v of the resulting LR tableau are easily 
recovered from these data. 

The combined equations (6.3) and (6.4) constitute a set of linear and homogeneous 
inequalities. We call this the LR (or tensor-product) basis. As described in [12], the Hilbert 
basis theorem guarantees that every solution can be expanded in terms of the elementary 
solutions of these inequalities. This is a key concept for the following (see [13] for an 
extensive discussion of these methods). A sum of two solutions translates into the product 
of the corresponding couplings, more precisely, to the stretched product (denoted by •) of 
the corresponding two LR tableaux. This is defined as follows. Denote the void boxes of 



a LR tableau by a 0, so that noj = J^iLj 1 A tableau is thus completely characterized 
by the data {riij} where now z > 0. Then, the tableau obtained by the stretched product 
of the tableaux {n^} and {n'^} is simply described by the numbers {riy + ?^-}, e.g., 
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(6.5) 



Let us now turn to the sit (2) case. The complete set of inequalities for sit(2) variables 
{Ai,nn,ni 2 } is simply 

Ai > ni2 tin > ni2 > (6.6) 

The first one expresses the fact that two boxes marked by a 1 cannot be in the same column 
while the other two are obvious. The other weights are fixed by the relation (ii = tin +ni2 
and v\ = Ai + nn — n\2- Any solution of these inequalities describes a coupling. By 
inspection, the elementary solutions of this set of inequalities are 



(Ai,n n ,n 12 ) = (1,0,1), (1,0,0), (0,1,0) 



(6.7) 



(For more complicated cases, we point out that powerful methods to find the elementary 
solutions are described in [10].) These correspond to the following LR tableaux, denoted 
respectively Ei : E 2 .>E 3 : 



Eo 



□ , 



E% 



(6.8) 



It is also manifest that there are no linear relations between these couplings. Any stretched 
product of these elementary tableaux is an allowed su(2) coupling. Because there are no 
relations between the elementary couplings, this decomposition is unique. We thus see 
that the description of the elementary couplings captures, in a rather economical way, the 
whole set of solutions of (6.6), that is, the whole set of sit(2) couplings. 



7. Reformulating the fusion rules in terms of linear inequalitites 

Consider now the affine-fusion extension of the reformulation of the sit(2) tensor prod- 
ucts in terms of linear inequalitites. The elementary couplings have a natural affine exten- 
sion, denoted by a hat, and their threshold levels are easily computed from the Kac- Walton 



formula.0 The result is: ko(Ei) = 1 for i = 1, 2,3. We observe that these values of ko are 
the same as the number of columns. Since the product of fusion elementary couplings is 
also a fusion and because this decomposition is unique, we can read off the threshold level 
of any coupling, hence of any LR tableau, simply from the number of its columns: 

ko = ^columns = Ai + n\\ (7.1) 

And since k is necessarily greater that ko, we have obtained the extra inequality: 

k>X 1 + n 11 (7.2) 

This together with (6.6) yield a set of inequalities describing completely the fusion rules. 
This is what we call a fusion basis, here the fusion basis of su(2). As in the finite case, the 
fusion couplings can be described in terms of elementary fusions. These correspond to the 
elementary solutions of the four inequalities, which are easily found to be 

(Mi,n n ,n 12 ) = (1,0,0,0), (1,1,0,1), (1,1,0,0), (1,0,1,0) (7.3) 

They correspond respectively to the coupling 



E Q : [1, 0] x [1, 0] D [1, 0] E 2 : [0, 1] x [1, 0] D [0, 1], 

E x : [0, 1] x [0, 1] D [1, 0] % : [1, 0] X [0, 1] D [0, 1]. 



(7.4) 



Any fusion has an unique decomposition in terms of these elementary couplings. For 
instance 



[3,2] x [1,4] D [1,4] 
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1 


1 





E 1 E 2 El : k = 5 (7.5) 



8. Constructing the fusion basis: Farkas' lemma 

For algebras other than su(2), the threshold level is not simply the number of columns. 
So the question is: how can we derive the fusion basis? The following strategy was devel- 
oped in [10]: 

2 Let us mention here that for the classical simple Lie algebras, the tableau methods for tensor 
products have been modified to implement the Kac- Walton formula for fusions - see [14] in the 
Sharp volume of the Canadian Journal of Physics. 



1) Write the LR inequalities; 

2) from these, find the tensor-product elementary couplings; 

3) from these, find fusion elementary couplings; 

4) from these, reconstruct the fusion basis. 

To go from step 2 to step 3, we need some tools; we describe below a method based on 
the outer automorphism group. Similarly to go from 3 to 4, we need a further ingredient: 
this is the Farkas' lemma. We discuss these techniques in turn. 

Let us start from the set of tensor-product elementary couplings {Ei, i G 1} for some 
set / fixed by the particular su(N) algebra under study. For each Ei, we calculate the 
threshold level ko(Ei) and this datum specifies the affine extension of Ei, denoted Ei. 
We have then a partial set of fusion elementary couplings with the set {Ei,i G /}. Our 
conjecture is that the missing fusion elementary couplings can all be generated by the action 
of the outer-automorphism group. For su(N), this group is simply {a n \ n = 0, . . . , N— 1}, 
with 

&[Ao, Ai, . . . , Ajv-i] = [Ajv-i, Ao, . . . , Aat_2] (8.1) 
The conjecture is based on the invariance relation 

a n X,a m p, A/2 v ' 

It amounts to supposing that the full set is contained in {AEi,i G I,WA}. Here AEi 
indicates a coupling of weights 

A{\, A; v} = {a n X, a m ji; a n+m P} , (8.3) 

n,m being arbitrary integers defined modulo N, if Ei has weights {A, /*;£}. The conjec- 
tured completeness requires the consideration of all possible pairs (n, m)JI 

Let us illustrate this with the su(2) case. Start with the elementary coupling E\ : 
(l)cg>(l) D (0), which, as already indicated, arises at level 1: ko(Ei) = 1. The corresponding 

3 Note that we do not suppose that the action of A on an elementary coupling will necessar- 
ily produce another elementary coupling. Indeed, the resulting coupling could be a product of 
elementary couplings. What is conjectured here is that all fusion elementary couplings can be 
generated in this way. 



fusion is thus [0, 1] x [0, 1] D [1, 0], denoted as E\. We now consider all possible actions of 
the outer-automorphims group on it. Since this group is of order 2, there are 4 possible 
choices for the doublet (n, m) : 

(a™, a m )e {(a, a), (1,1), (1, a), (a, 1)} (8.4) 

with a[Ao,Ai] = [Ai,Ao]. This generates the set of four elementary couplings found pre- 
viously, in the respective order Eq, Ei, E 2 , E 3 . Thus, from one tensor-product elementary 
coupling, all four fusion elementary couplings are deduced. 

We now turn to Farkas' lemma. For its presentation, it is convenient to use an expo- 
nential description of the couplings, that is, 

(fc^ny) - d k L^N^ (8.5) 

d, Li, Nij being dummy variables. For instance E\ is represented by dL\N\i- If we 
collectively describe a coupling by the complete set of variables {x^}, we have 

{ Xi } - (8.6) 

A particular coupling is thus described by a given product rii^T* with fixed Xi. Its 
decomposition in terms of elementary couplings take the form Yli E z ai . Now, since Ei can 
be decomposed in terms of the X, as 

3 

it means that reading off a particular coupling means that we are interested in a specific 
choice set of positive integers {xi} fixed by 

^ ai€ij = Xj (8.8) 

i 

in terms of some positive integers a^. We are thus looking for the existence conditions 
for such a coupling, i.e., the underlying set of linear and homogeneous inequalities. This 
is exactly what the Farkas' lemma [15] gives us: given the knowledge of the elementary 
couplings, it allows us to recover the underlying set of inequalities. For tensor products, 
this is of no interest since we know the corresponding set of inequalities and our elementary 



couplings have been extracted from them. But the situation is quite different in the fusion 
case, where the fusion basis is unknown. 

For our application we need the following modification of the lemma, proved in [10]: 

Lemma: Let A be an r x m integer matrix and let €j (j = 1, . . . , n) be a set of fundamental 
solutions to 

Ax > 0, xe N m . (8.9) 

Let V be the m x n matrix with entries Vij = (for i = 1, . . . , m, j = 1, . . . , n), i.e., 
the columns of V are a set of fundamental solutions to (8.9). Let e w (w = 1, ...,£) be a 
fundamental system of solutions of u T V > 0, (not necessarily positive) u G Z m , and let 
E be the t x m matrix with entries E wi = (e w )i, i.e., the rows of E are the fundamental 
solutions e w (w = 1, . . . , £). Then the solution set of the system 

Ex > 0, xeN m (8.10) 

is the same as the solution set of (8.9). 

To link the lemma to the situation presented above, we note that the entries of 
the matrix V are given here by the numbers €ji appearing in (8.7). The relation V a = 
x describes a generic coupling Yli &i ai ? an d our goal is to find the defining system of 
inequalities for x that underly the existence of this coupling. 

Take a simple example, the su(2) case. The elementary couplings and the correspond- 
ing vectors are 

E :d: e = (1,0,0,0), E 2 : dL x : e 2 = (1,1,0,0) 

E 1 :dL 1 N 12 : e 1 = (1, 1, 0, 1) , E 3 : dNu : e 3 = (1, 0, 1, 0) 

From the vectors €i with components eij, we form the matrix V with entries Vij = e^: 



V = 



/I 1 1 1\ 

110 

1 

Vo i o o/ 



(8.12) 



With a and a; denoting the column matrices of entries ai and X{ respectively, we have the 
matrix equation 

Va = x (8.13) 



Again, this equation describes a general fusion coupling Yl i E t ai . We now want to unravel 
the underlying system of inequalities. For this, we consider the fundamental solutions of 

u T V>0 (8.14) 

where u is the vector of entries Ui. These inequalities read 

uq > , uq + u\ + 113 > , uq + ui > , uq + U2 > (8.15) 

In this simple case, the elementary couplings can be found by inspection and these are: 

e = (1,-1,-1,0), e 1 = (0,0,0,1), e 2 = (0,1,0,-1), e 3 = (0,0,1,0) (8.16) 

Finally, we consider the conditions ej x > 0, with (xo, xi, X2, X3) = (k, Ai, nn, 77-12). They 
read, in order, 

k > Ai + nn, ^12 > 0, Ai > ni2, nn > (8.17) 

The last three conditions define the LR basis. The first one is the additional fusion con- 
straint. Together, they form the su(2) fusion basis. 

9. Constructing the fusion basis: polytope techniques 

In the previous section, the Farkas' lemma has been used to construct the fusion basis 
out of the set of fusion elementary couplings. There are alternative approaches, however. 
Another one is based on the reinterpretation of the fusion-rule computations in terms of 
counting points inside a polytope. A polytope can be described by its vertices or its facets. 
In our context, the vertices are represented by the fusion elementary couplings and the 
facets are the inequalities for which the elementary couplings are the elementary solutions. 
The reconstruction of the facets of a polytope from its vertices is thus another way to 
generate the fusion basis. This method is described in [16]. 

In the special case of su(N), Berenstein-Zelevinsky triangles can also be used to derive 
the polytope description of a fusion basis [17], by considering so-called virtual triangles. 
Multiple sum formulas can then be written for fusion coefficients of various types. However, 
as are all methods to date, this one is difficult to extend to higher rank. Assigning a 
threshold level to a Berenstein-Zelevinsky triangle becomes very rapidly more difficult 
with increasing rank (see [18]). 



10. Conclusion 



The fusion bases have been constructed for su(3), su(4) and sp(4). Note that for 
algebras other that su(N), we replace the LR basis by the Berenstein and Zelevinsky 
basis [19]. This leads to explicit expressions for the threshold levels, hence for the fusion 
coefficients.il 

We stress that the reformulation of the problem of computing fusion rules in terms of 
a fusion basis solves, in principle, the quest for a combinatorial method, since it reduces a 
fusion computation to solving inequalities. But it is probably not the optimal solution to 
the quest for an efficient combinatorial description. 

The main open problem concerning fusion bases is to find a fundamental and Lie alge- 
braic way of deriving the basis, analogous in spirit to the Berenstein-Zelevinsky conjectures 
for generic Lie algebras in [19]. 

The methods described in section 8, involving Farkas' lemma, are general and powerful, 
but they may not be Lie algebraic enough. Perhaps one should step back and look at a 
first principles description of the tensor product couplings, and its adaptation to fusion. 
This was done in [20]. Three-point functions were calculated that can be regarded as 
generating functions for tensor product couplings, and a very simple method was found 
for adapting the results to fusion couplings. In principle, the procedure works for any 
semi-simple Lie algebra. Unfortunately, these more Lie algebraic methods are inevitably 
more involved. It gives more information (such as operator product coefficients instead of 
just fusion coefficients) , but it is not clear that it can be implemented effectively on higher 
rank algebras. 



Related methods for obtaining the threshold level are presented in [18]. 
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